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(54) Method and device for assessing the stability of an electric power transmission network 



(57) In a method and devices for assessing the sta- 
bility of an electric power network, for each of a plurality 
of measurement locations with associated voltage and 
current measurements, a normalized power margin val- 
ue dSn is computed, which is defined as a result of di- 
viding an associated power margin value by an associ- 
ated maximum allowable power flow. In a preferred em- 
bodiment of the invention, the normalised power margin 
value dSn is computed from an apparent load imped- 
ance 2a and a Thevenin impedance Zt, that are per- 



ceived from the measurement location, as 



— esr- 



Normalised power margin values from several 
measurement locations in the network therefore show 
power margins relative to a local maximum capacity, 
which makes them comparable to one another. 
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Description 

Field of th Inv ntlon 

[0001] The invention relates to large-scale electric 
power transmission networks, and, more particularly, to 
a method and a computer program product for assess- 
ing the stability of an electric power network, and a pow- 
r network stability assessment device according to the 
preamble of claims 1 , 7, 8 and 11 , respectively. 

Background of the Invention 

[0002] Electric power transmission and distribution 
systems or networks comprise high-voltage tie lines for 
connecting geographically separated regions, and sub- 
stations for transforming voltages and for switching con- 
nections between lines. Power generation and load flow 
in a network with several substations is controlled by a 
central energy management system. An important issue 
in the control of a power generation and load flow is to 
keep the network stable, i.e. to avoid voltage collapses 
and swings. A method for assessing network stability, 
based on voltage margins, is described in the paper 
"Use of local measurements to Estimate Voltage-Stabil- 
ity Margin", K. Vu et al., Power Industry Computer Ap- 
plications (PICA) May 12-16, 1997, IEEE, and in "Volt- 
age instability predictor (VIP) and its applications", K. 
Vu et al., Power Systems Computation Conference 
(PSCC) June 1999. Both articles are herewith incorpo- 
rated by reference. These articles describe a "Voltage 
Instability Predictor" (VIP) which measures currents and 
voltages locally in order to infer a proximity to voltage 
collapse. The concept of the VIP is shown in Figure 1 . 
One part of an electric power system is treated as a pow- 
er source, another part as a load. The power source is 
r presented by its TheVenin equivalent 1 with a source 
voltage ft and a Thevenin or source impedance ?t. The 
load is represented by an apparent load impedance 2*a. 
Both the Thevenin impedance 2t and the apparent load 
impedance ?a are estimated from the current and volt- 
age measurements by a VIP device 2. The relation of 
these impedances, expressed by a stability margin or 
power margin, indicates how close the power system or 
n twork is to collapsing. 

[0003] The power margin indicates, in MVA (Mega 
Volt Ampere), how much power can be drawn from a 
substation or transmitted through a tie line until the volt- 
age collapses. If there is a plurality of VIP devices 2, it 
is necessary to process their associated power margin 
values in a manner that gives a meaningful combined 
representation of a state of the network. Comparing the 
power margin values themselves is not meaningful. 
That is, a power margin value of 10 MVA has another 
meaning for e.g. a transmission line with a rated thermal 
capacity of 5 kA (5000 Amperes) than in a line with a 
rated capacity of 1 kA. 



Description f th Inv ntion 

[0004] It is therefor an object of th inv ntion to cre- 
ate a method and a computer program product for as- 

s sessing the stability of an electric power network, and a 
power network stability assessment device of the type 
mentioned initially, which allow to determine local net- 
work stability margin information at a plurality of meas- 
urement locations of the network, such that said infor- 

10 matlon from different measurement locations is compa- 
rable in a meaningful manner. 

[0005] These objects are achieved by a method and 
a computer program product for assessing the stability 
of an electric power network, and a power network sta- 
15 bility assessment device according to the claims 1 , 7, 8 
and 11, respectively. 

[0006] In the inventive method, computer program 
product and device, for each of a plurality of measure- 
ment locations with associated voltage and current 

20 measurements, a normalized power margin value is 
computed, which is defined as a result of dividing an as- 
sociated power margin value by an associated maxi- 
mum allowable power flow. Both the power margin value 
stability margin and the maximum allowable power flow 

25 depend on the measured values. 

[0007] Normalized power margin values from several 
measurement locations in the network therefore show 
power margins relative to a local maximum capacity, 
which makes them comparable to one another. This has 

30 the advantage that, in a critical situation, when the net- 
work approaches instability, a corresponding cause can 
be localized from locations of measurements that show 
a low normalized power margin. A further advantage is 
that the local maximum capacity is a function of an ac- 

35 tual state of the network and varies according to the 
measurements. It therefore reflects the state of the net- 
work better than a rated maximum value, which remains 
a constant design value unrelated to reality. 
[0008] In a preferred embodiment of the invention, 

40 normalized power margin values from the plurality of 
measurement locations in the network are represented 
together as a bar graph in a graphic display. 
[0009] In a preferred embodiment of the invention, the 
power network stability assessment device is a power 

45 network stability comparison device that is configured 
to receive information from a plurality of measurement 
devices that measure voltage and current values at a 
plurality of measurement locations in the electric power 
distribution network. The power network stability com- 

50 parison device comprises means to determine, for each 
of the plurality of measurement locations, a normalized 
power margin value that is defined as a result of dividing 
an associated power margin value by an associated 
maximum allowable power flow. 

55 [0010] In another preferred embodiment of the inven- 
tion, the power network stability assessment device is 
a power network stability estimating device that com- 
prises means for measuring voltage and curr nt values 
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at a location in the electric power distribution network 
and means for estimating an apparent load impedance 
magnitud Za and an estimat d Thevenin impedance 
magnitude Zt associated with th measurement loca- 
tion. It comprises means for computing a normalised 5 
power margin value dSn that is essentially equal to 



dSn-fS=»? 
Ua+zt; 



[0011] Further preferred embodiments are evident 
from the dependent patent claims. 

Brief Description of the Drawings 

[001 2] The subject matter of the invention will be ex- 
plained in more detail in the following text with reference 
to preferred exemplary embodiments which are illustrat- 
ed in the attached drawings, in which: 

Figure 1 shows schematically a conceptual 

structure for assessing network stabil- 
ity, according to the state of the art; 

Figure 2 is a diagram showing voltages versus 

current at a selected point in an elec- 
tric power network; 

Figure 3 schematically shows a power network 

with measurement and data process- 
ing devices according to the invention; 

Figures 4 to 6 show a display of normalized power 
margins according to a first variant of 
the invention; and 

Figure 7 shows a display of normalized power 

margins according to a second variant 
of the invention. 

[0013] The reference symbols used in the drawings, 
and their meanings, are listed in summary form in the 
list of reference symbols. In principle, identical parts are 
provided with the same reference symbols in the figures. 

Detailed Description of Preferred Embodiments 

[0014] Figure 1 schematically shows a conceptual 
structure for assessing network stability, in which one 
part of an electric power system is treated as a power 
source, another part as a load. The power source is rep- 
resented by its Thevenin equivalent 1 with a Thevenin 
or source impedance 2t. The load is represented by an 
apparent load impedance 2a. A VIP device 2 deter- 
mines phasor data with a phasor measurement unit re- 
siding, for example, at a feeder at the bay level of sub- 
stations or at branching points along transmission lines. 
A voltage phasor represents, for example, a voltag of 
the feeder or line, whil a current phasor represents cur- 



rent flowing through th f eder or lin . 
[001 5] The phasor data represents a phasor and may 
b a polar number, th absolute value of which corre- 
sponds to either the real magnitude or the RMS valu 
of a quantity, and the phase argument to the phase an- 
gle at zero time. Alternatively, the phasor may be a com- 
plex number having real and imaginary parts or the 
phasor may use rectangular or exponential notation. 
Phasors may be used to represent quantities such as 
the voltage, current, power or energy associated with a 
phase conductor or an electronic circuit. By contrast, 
conventional sensing devices used in power networks 
generally measure only scalar, average representa- 
tions, such as the RMS value of a voltage, current etc. 
[0016] In some VIP applications, the phasor data is 
collected from phasor measurement units that are dis- 
tributed over a large geographic area, i.e. over tens to 
hundreds of kilometers. For applications in which the 
phasor data from these disparate sources are analysed 
in conjunction, they must refer to a common phase ref- 
erence. In other words, the different phasor measure- 
ment units must have local clocks that are synchronised 
to within a given precision. Such a synchronisation of 
the phasor measurement units is preferably achieved 
with a known time distribution system, for example the 
global positioning (GPS) system. In atypical implemen- 
tation, the phasor data 9 is determined at least every 
200 or every 100 or preferably every 40 milliseconds, 
with a temporal resolution of preferably less than 1 mil- 
lisecond. In a preferred embodiment of the invention, the 
temporal resolution is less than 1 0 microseconds, which 
corresponds to a phase error of 0.2 degrees. Each 
measurement is associated with a time stamp derived 
from the synchronised local clock. The phasor data 
therefore comprises time stamp data. 
[0017] The VIP assesses the stability of the electric 
power transmission network by determining a stability 
margin value pertinent to specific entities and/or a com- 
bination of entities within the network. A measure of sta- 
bility is expressed in terms of impedances or, as a volt- 
age margin, in terms of voltages. In the example de- 
scribed in the background of the invention, a power mar- 
gin is used as a stability margin. A power margin is a 
more intuitive representation than a voltage or imped- 
ance margin. A local power margin represents for ex- 
ample the amount of power that may be transmitted 
through a given transmission line of the network. A glo- 
bal power margin combines phasor data collected from 
a plurality of phasor measurement units. 
[001 8] Figure 2 shows a plot of voltage versus current 
at a given point in a power network, where a VIP device 
is located. This may be at a feeder of a tie line or a load. 
When the load gets stronger, the current increases and 
the voltage drops. The actual future voltage versus cur- 
rent curve (V/l curve) is represented by the dashed line 
4. The actual future behavior is influenced by the entire 
network, but is estimated from local measurements by 
th VIP device. The VIP computes provides an estimate 
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of this curv , shown by th drawn out line 3. In particular, 
th slope of the curve corresponds to the magnitude ifti 
ofth Thev ninimpedanc Zt computed by the VIP d - 
vice. In the r maining text, said magnitud shall be d - 
noted as Zt, and the magnitude of the apparent load im- 
pedance ta as Za. Note that in the following computa- 
tion of impedances, the phasors as well as the imped- 
ances are vector quantities. The computation of voltage 
margins, shown further on and in Figure 2, involves cor- 
responding scalar values, i.e. magnitudes of impedanc- 
s, voltages and currents. 
[0019] The apparent load impedance 2a is deter- 
mined from a first measurement comprising a first volt- 
age phasor \?1 and a first current phasor Tl , a first load 
impedance is computed as 



fa1 = 



\T1 

TP 



[0020] The Thevenin impedance 2t is estimated from 
th first measurement and at least one second meas- 
urement comprising a second voltage phasor V*2 and a 
second current phasor T2, 

* f V*2-v*1 

[0021] Each of the measurements represents an op- 
erating point on the actual V/l curve 4. 
[0022] In order to improve an accuracy of the esti- 
mate, statistical and/or heuristical methods as shown, 
for example, in US patent 5,631,569, are used. The 
magnitude of the Thevenin impedance Zt is estimated 
by taking advantage of natural fluctuations in the power 
network, which cause measured voltages and currents 
to change. A RMS (root mean square) voltage and an 
RMS current are measured cyclically. For each cycle, 
an apparent load impedance is computed. The Theven- 
in impedance is only computed if two cycles with dispa- 
rate load impedances are identified. In a preferred em- 
bodiment of US patent 5,631,569, a 10% difference in 
load impedance is considered useful for calculating 
Thevenin impedance. Successive values of the Theven- 
in impedance are stored and statistical data is generat- 
ed and maintained for the successive values. When a 
standard deviation falls below one sigma, the mean of 
the values is displayed. 

[0023] According to the invention, given a first oper- 
ating point 5 with voltage magnitude Vp, current magni- 
tude Ip and estimated Thevenin impedance magnitude 
Zt, any other operating point with voltage V and current 
I on the estimated V/l curve 3 satisfies the equation 

(V-Vp)=-Zt(Mp). 



40 



45 



50 



[0024] The maximum power corresponds to an sti- 
mated maximum power operating point 6 on th esti- 
mated V/l curve 3 for which th product of curr nt and 
voltag is maximal. This product is qual to 



Vxl=(Vp-ZMp)l-ZW 



and the maximum of the product is equal to 



10 



15 



20 



25 



30 



35 



Smax= <VP+Z"p)* 
4Zt 



[0025] This maximal power is the power transmitted 
at the maximum power operating point 6. It gives the 
maximum power, in MVA, that may be transmitted 
through said feeder or line under current network con- 
ditions. 

[0026] The power margin, for a present operating 
point 5, is the difference between the maximum power 
and the power being currently transmitted. Let the 
present time current and voltage magnitudes be I and 
V, respectively. Then the present power margin or 
present stability margin dS corresponding to said feeder 
Is 



dS = 



(V+ZMr 
4Zt 



•Vxl 



which is equal to 



dS = 



4-Zt 



[0027] In summary, the stability margin computation 
method gives, based on online current and voltage 
phasor measurements, a continuously adapted esti- 
mate of how much more power may be transferred 
through a tie line or how much more power may be 
drawn by a substation before the network collapses. 
This is particularly advantageous in situations where the 
network state slowly moves towards instability, without 
tell-tale disruptive events that would indicate a critical 
situation. 

[0028] The normalised power margin value dSn ac- 
cording to the invention is 



dSn: 



dS _ (V-ZMr 



Smax 



(V+ZMr 



[0029] Dividing the numerator and denominator by I, 
and given that Za = V / 1, the following result is obtained 
for the normalized power margin 
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[0030] This normalized power margin dSn will vary s 
from 1 00% at no load to 0% at voltage collapse. Regard- 
less of the VIP location, all normalized power margins 
can be easily compared to determine which locations 
are closer to voltage collapse. 

[0031] Under normal conditions, the apparent load 10 
impedance magnitude Za is much greater than the 
Thevenin impedance magnitude Zt. As load increases, 
apparent load impedance Za correspondingly decreas- 
es, while the Thevenin impedance Zt may increase. 
When the impedances become equal, maximal power 15 
is delivered to the load and the system is on the verge 
of collapse. 

[0032] If Za decreases below Zt, the system becomes 
voltage unstable and power delivered to the load actu- 
ally decreases. Under this extreme condition, the term 20 
Za - Zt will be a negative number, but the above equation 
for dSn will still evaluate to a positive number. Therefore, 
In a preferred embodiment of the invention, in order to 
indicate a system condition clearly, the equation is mod- 
ified as follows so that the normalized power margin is 2s 
negative in the unstable region: 

dSn- = ^(Za-Zt)-(g^) 2 , 

where s/on(Za-Zt) is the signum function. Its value is 1 
if the bracketed term is larger than or equal to zero, and 
-1 if it is smaller than zero. 

[0033] In summary, the inventive method for assess- 35 
ing the stability of an electric power transmission or dis- 
tribution network, where said network comprises a plu- 
rality of substations, buses and lines, comprises the step 
of 

40 

a) measuring voltage and cu rrent values V,l at a plu- 
rality of measurement locations in the electric power 
network 31 ,32,33, and 

b) determining, for each of the plurality of measure- 
ment locations, a normalized power margin value 45 
dSn that is defined as a result of dividing an asso- 
ciated power margin value dS by an associated 
maximum allowable power flow Smax. 

[0034] In different preferred variants of the invention, so 
the normalised power margin value dSn is computed 

■ from measured voltages and currents and/or 

■ from intermediate variables such as the Thevenin 
impedance 2t and apparent load impedance 2a or 55 

■ from their respectrv magnitudes Zt.Za or 

■ from th power margin value dS and the maximum 
allowable power flow Smax. 



[0035] Thes variants are mathematically essentially 
equal, as are results obtained through the drff rent var- 
iants. As such, they all fall und r the scope of the inven- 
tion. 

[0036] Figure 3 schematically shows a power network 
with measurement and data processing devices accord- 
ing to the invention. The power network comprises sub- 
stations 31 connected by tie lines 32 and comprising 
busbars 33 and measurement devices 34. The meas- 
urement devices 34 are VIP devices 2 and/or phasor 
measurement units. They are configured to communi- 
cate through communication links 35 with a central unit 
36 or power network stability comparison device 36. The 
central unit 36 comprises a display unit 37 or is config- 
ured to communicate with a display unit 37. 
[0037] In a first preferred embodiment of the inven- 
tion, measurement devices 34 transmit phasor informa- 
tion, that is, measured and time stamped voltage and 
current information through the communication links 35 
to the central unit 36. The central unit 36 thus comprises 
means for determining, for connected measurement de- 
vices 34, the normalised power margin value dSn from 
their measurement values. For this first embodiment, it 
is sufficient that the measurement devices 34 are 
phasor measurement units. 

[0038] In a second preferred embodiment of the in- 
vention, measurement devices 34 determine and trans- 
mit time stamped information about values of the appar- 
ent load impedance 2a and the estimated Thevenin im- 
pedance 2t, or about their respective magnitudes. Ac- 
cordingly, the central unit 36 comprises means for re- 
ceiving, from connected measurement devices 34, said 
impedance information through the communication 
links 35 and means for computing the normalised power 
margin value dSn from the impedance information. 
[0039] In a third preferred embodiment of the inven- 
tion, measurement devices 34 determine and transmit 
time stamped information about values of the normal- 
ised power margin value dSn. Accordingly, the central 
unit 36 comprises means for receiving the normalised 
power margin value dSn through the communication link 
35. 

[0040] In a fourth preferred embodiment of the inven- 
tion, measurement devices 34 determine and transmit 
time stamped values of phasors and of the Thevenin im- 
pedance Zt or its magnitude Zt, from which the central 
unit 36 computes the normalised power margin value 
dSn. 

[0041] In a preferred embodiment of the invention, the 
normalized power margin values from the plurality of 
measurement locations and corresponding locations in 
the network are represented together in a graphic dis- 
play, in particular as a bar graph. 
[0042] Figures 4, 5 and 6 show a display of normal- 
ized power margins dSn according to a first variant of 
the invention. The normalized power margins dSn from 
a plurality of numbered measurement locations ar 
shown along a horizontal axis. Their respective values 
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are drawn along a vertical axis from 0 to 100%, with val- 
ues indicated numerically and graphically by a height of 
a corresponding bar. Figure 4 shows all normalized 
power margins dSn being fairly larg , indicating that 
there is no imminent danger to network stability. In figure 
5, stability margins throughout the network are lower 
than in figure 4, indicating that network stability is in dan- 
ger. The measurement locations that lead to the lowest 
values of the normalized power margins dSn are of spe- 
cial interest: they indicate what part (in the topologicaJ 
or geographical sense) of the network is in a most critical 
state. This information helps an operator to plan and pri- 
oritize countermeasures such as load shedding or di- 
verting load flow. In figure 5, the power network is on the 
verge of collapse. 

[0043] Figure 7 shows a display of normalized power 
margins dSn according to a second variant of the inven- 
tion. Vertical bars are associated with measurement lo- 
cations, but here a distance from the horizontal 100% 
line indicates the local power margin. In other words, the 
height of a bar is proportional to the amount of power 
corresponding to the current operating point 5, divided 
by the maximum allowable power flow Smax. The net- 
work state corresponding to Figure 7 is the same as for 
Figure 6. 

[0044] A computer program product according to the 
invention comprises a computer readable medium, hav- 
ing thereon computer program code means to make, 
when said program is loaded, the computer execute the 
method according to one of the preceding claims. 

List of designations 



[0045] 

1 

2 

3 
4 
5 
6 

31 
32 
33 
34 
35 
36 

37 
V*1 
V*2 

T2 
T2 
Vp 



Thevenin equivalent 
voltage instability pre- 
dictor, VIP device 
estimated V/l curve 
actual future V/l curve 
operating point 
maximum power oper- 
ating point 
substation 
tie line 
busbar 

measurement device 
communication link 
central unit, compari- 
son device 
display unit 
first voltage phasor 
second voltage 
phasor 

first current phasor 
second current phasor 
voltag magnitude at 
operating point 



w 



15 



20 



25 



30 



35 



40 



45 



IP 

2*a 

Za 

ft 
Zt 

dS 
dSn 

Smax 



Claims 



current magnitud at 
operating point 
apparent load imp d- 
anc 

apparent load Imped- 
ance magnitude 
Thevenin impedance 
Thevenin impedance 
magnitude 
power margin value 
normalised power 
margin value 
maximum allowable 
power flow 



2. 



3. 



50 



55 



Method for assessing the stability of an electric 
power network, where said network comprises a 
plurality of substations, buses and lines, and 
where the method comprises the step of 

a) measuring voltage and current values V.l at 
a plurality of measurement locations in the elec- 
tric power network (31,32,33), 
characterised in that the method further com- 
prises the step of 

b) determining, for each of the plurality of meas- 
urement locations, a normalized power margin 
value dSn that is defined as a result of dividing 
an associated power margin value dS by an as- 
sociated maximum allowable powerflow Smax, 
where both the power margin value stability 
margin dS and the maximum allowable power 
flow Smax depend on the measured values V,l. 

Method according to claim 1, characterised in 
that, for each of the plurality of measurement loca- 
tions, an apparent load impedance Za and an esti- 
mated Thevenin impedance Zt are estimated and 
that the normalised power margin value dSn is com- 
puted as being essentially equal to 



— ear- 



Method according to claim 1 , characterised in that 
the maximum allowable powerflow Smax is defined 



4-Zt 



where V is a measured current, I is a measured volt- 
age and Zt is an estimat d Thevenin impedance of 
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th network. 

4. Method according to claim 1 , characterised in that 
the power margin value dS is defined as 

5 

ds= (V-ZM) 2 
42t 

where V is a measured current, I is a measured volt- 10 
age and Zt is an estimated Thevenin impedance of 
the network. 

5. Method according to claim 1 , characterised in that 
normalized power margin values from the plurality is 
of measurement locations in the network are repre- 
sented together in a graphic display. 

6. Method according to claim 5, characterised in that 
normalized power margin values from the plurality 20 
of measurement locations in the network are dis- 
played together in a bar graph. 



1 1 . Power network stability assessment device, where- 
in the stability assessment device is a stability esti- 
mating device that comprises means for measuring 
voltage and current values at a location in an elec- 
tric power network (31,32,33), and means for esti- 
mating an apparent load impedance Za and an es- 
timated Thevenin impedance Zt associated with the 
measurement location, characterised in 
that it comprises means for computing a normalised 
power margin value dSn that is essentially equal to 



7. A computer program product comprising computer 
program code means to make, when said program 25 
is loaded, the computer execute the method ac- 
cording to one of the preceding claims. 



8. Power network stability assessment device, where- 
in the assessment device is a comparison device 30 
(36) that is configured to receive information from a 
plurality of measurement devices (34) that measure 
voltage and current values at a plurality of meas- 
urement locations in an electric power network 
(31,32,33), characterized in that the power net- 35 
work stability comparison device (36) comprises 
means to determine, for each of the plurality of 
measurement locations, a normalized power mar- 
gin value dSn that is defined as a result of dividing 
an associated power margin value dS by an asso- *o 
dated maximum allowable power flow Smax. 



9. Device according to claim 8, characterised in that 
the means to determine, for each of the plurality of 
measurement locations, the normalised power mar- 45 
gin value dSn, comprise means for receiving the 
normalised power margin value dSn through a com- 
munication link (35). 



1 0. Device according to claim 8, characterised in that so 
the means to determine, for each of the plurality of 
measurement locations, the normalised power mar- 
gin value dSn, comprises means for receiving infor- 
mation about an estimated apparent load imped- 
ance Za and an estimated Thevenin impedance Zt 55 
through a communication link (35) and means for 
computing the normalised power margin valu dSn 
that is essentially equal to 
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